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We investigate the uncertainties of the heavy-quark parton distribution functions in the variable 
flavor number scheme. Because the charm- and bottom-quark parton distribution functions (PDFs) 
. . . are constructed predominantly from the gluon PDF, it is a common practice to assume that the 

• heavy-quark and gluon uncertainties are the same. We show that this approximation is a reasonable 

' first guess, but it is better for bottom quarks than charm quarks. We calculate the PDF uncertainty 

^— , for t-channel single-top-quark production using the Hessian matrix method, and predict a cross 

C^l . section of 2.12toll Pb at run II of the Tevatron. 

> ; 

Q , As a new run of the Fermilab Tevatron begins, there is a considerable interest in measuring fundamental 
parameters of the Standard Model and in looking for new particles. For some of these measurements the 
dominant uncertainty will come from parton distribution functions (PDFs). In particular, measurements of the 
Cabibbo-Kobayashi-Maskwawa (CKM) matrix element Vtb depend on an accurate prediction of the cross section 
' for single-top-quark production and will be limited by uncertainties in the &-quark PDF ||^ . In a large tan f3 
I , model of supersymmetry, Higgs production from a bb initial state is the dominant production mechanism . 
^ ■ In order to understand uncertainties of the heavy-quark parton distribution functions, it is necessary to 
00 remember where they arise. In the conventional QCD framework, heavy quarks are produced perturbatively 
from gluons splitting into heavy-quark pairs at energies above the heavy-quark mass threshold (we do not 
CO I consider "intrinsic" heavy quarks Q here). The heavy-quark PDF Q(x,fi^) is a formal object, which can be 
used reliably when the typical momentum scale /i of the hard scattering is much larger than the mass of the heavy 
quark mq. This PDF resums nearly collinear singularities a" In" (/i^/mQ)/n! associated with propagation of the 
initial-state heavy quarks. Unlike the light-quark PDFs, the distributions for charm (c) and bottom quarks (b) 
are not independent functions that are determined by fitting to the hadronic data; instead, they are generated 
_ , from the PDFs of light partons in the process of Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution. 
II ' The relationship between the PDFs for heavy and light partons is especially simple near the mass threshold 
Oh, of heavy quarks (/i ~ mg), where the DGLAP equations can be solved in the leading-logarithm approximation. 
^ ' Given the initial condition Q{x, y?) = at /x = mg, and neglecting gluon bremsstrahlung off heavy-quark lines 
and the scale dependence of the gluon PDF and a^, the leading-order solution can be found as [5-7] 



X 




g(x,M^) = ^ln ^Pg, ^ (1) 



where the function PQg{y) = [y^ -I- (1 — y)^] /2 describes the splitting of a gluon into a heavy-quark pair. 
Since g{x) grows as (n « 1.4-1.5), at sufficiently small x Eq. ([^) may be rewritten as 

Til " f-- "<>M^) p) 

g{x,y^) as[iL'^) \n n+1 n + 2J \niQ J \mQ J 

The approximation (|^) works remarkably well over a wide range of x and /i. In Fig. ^ we present an updated 
version of Fig. 5 of Ref. ||l|. The ratios Q{x, /i^)/ g{x, /x^) x 27r/as(/^^) are shown as functions of /i for various fixed 
values of x, using the CTEQ5M1 parton distribution functions S. The dependence on ln(/i/r7iQ) is approxi- 
mately linear indicating that, even at next-to-leading order (NLO), Q(x,y^) oc [as(/i^)/27r] hi{y'^ /niQ)g{x, y"^). 
Further, the constant of proportionality saturates at ^ 0.5 for charm quarks when x ^ 0.1, and for bottom 
quarks when x < 0.05. 
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FIG. 1: The ratio of the c and 6 distribution functions to the gluon distribution function, times 27r/aa(At'^), versus 
ln(/i/mQ) for various fixed values of x. The curves are approximately linear, while the slope of the curves saturates at 
about 0.5 at small x, in agreement with the approximation of Eq. (^). 

If the heavy-quarks PDFs are directly proportional to the gluon PDF, it is reasonable to expect that their 
uncertainties are approximately the same: 



Q{x,iJ?) 



(3) 



where kq « 1. We test this relation by estimating the PDF uncertainties with the Hessian matrix method 
proposed J. Pumplin, et al. This paper contains 16 pairs of PDF sets, where each pair corresponds to 

varying one independent parameter Zj in the PDF fit such that the of the fit changes by — (5)^ = 25. We 
define the maximum positive and negative errors on an operator O by P, 0] 
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^(max[O(z,f + <)-O(z0),O(z0-t)-O(z0),0]) , 



2 

^(max[O(z0)-O(z0 + 0,Ok°)-O(^."-i),0]) , 



(4) 
(5) 



where the "tolerance" T is a scaling parameter that determines the overall range of allowed variation of |^]. 
In this study we use T = 10. 

Figure ^ shows the dependence of the ratios Kc and K}, on x and /i. We see that these ratios are quite close 
to unity for x < 0.05, but climb to about 2 at larger x. The approximation (^) holds better for h quarks: is 
less than 1.4 at a; < 0.1 and the whole range of /i, while Kc can reach up to 1.75 in this region. The values kq 
are close to 1 as a; —> 0, which supports the validity of Eq. (^) in the small- a; region. 

Figure ^ shows the relative uncertainties 5c{x, fj?)/c{x, y?) and 5b{x, iJ,'^)/b{x, //^) as functions of x at various 
values of /x. Both uncertainties have minima around x w 0.01, where the PDFs are best constrained by the 
existing data. At small x, the uncertainties grow because the small-x region is covered only by DIS experiments, 
which do not constrain well the gluon PDF. At x > 0.1, the heavy-quark PDFs become negligible compared to 
the valence quark PDFs, so that c(x,/i^) and b{x,y?) are practically unconstrained at a; ^ 0.3. 

Single-top-quark production via t-channel W^-exchange probes the 6-quark PDF at /i ^ mt ~ 175 GeV. The 
range of x probed at the run II of the Tevatron for accepted events is 0.06-0.5, with the bulk of the cross section 
coming from the region around x ~ 175/2000 ~ 0.09. Using this value of x and Fig. we can estimate the 
uncertainty of the total cross section at run II to be ^ ±15%. We have explicitly re-calculated the NLO cross 
section in Ref. and its uncertainty with the method described above and T = 10. The result is 2.12^Q'2g pb, 
or ti4%, in excellent agreement with Fig. |3|d. For this choice of T, and using CTEQ5M1 PDFs, the PDF 
uncertainties will be larger than all other theoretical or experimental errors once 2 fb^^ of integrated luminosity 
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FIG. 2: The ratios hq = {SQ /Q) /{5g/ g) for c and b quarks as a function of x and fi. 
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FIG. 3: The uncertainties of the c and b distribution functions as a function of x for various values of scale /i. 

is accumulated. One positive note is that Fig. ||b predicts and uncertainty of around ±7% at the LHC, which 
will be comparable to experimental systematics. 

The approximation that the heavy-quark uncertainties are the same as the gluon uncertainty is a good 
guess for smaller x and /i. We calculate these uncertainties using a "tolerance" T = 10 in a Hessian matrix 
method We show that the heavy-quark uncertainties strongly depend on x and relatively weakly on the 
scale fi. We calculate the PDF uncertainty for t-channel single-top-quark production, and predict a cross section 



of 2.12: 



-0.29 



pb at run II of the Tevatron. 
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